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ABSTRACT 

We use a set of twelve high-resolution N-body/hydrodynamical simulations in the ACDM cosmology to inves- 
tigate the origin and formation rate of fossil groups (FGs), which are X-ray bright galaxy groups dominated 
by a large elliptical galaxy, with the second brightest galaxy being at least two magnitudes fainter. The simu- 
lations invoke star formation, chemical evolution with non-instantaneous recycling, metal dependent radiative 
cooling, strong star burst driven galactic super winds, effects of a meta-galactic UV field and full stellar pop- 
ulation synthesis. We find an interesting correlation between the magnitude gap between the first and second 
brightest galaxy and the formation time of the group. It is found that FGs have assembled half of their final 
dark matter mass already at z > 1, and subsequently typically grow by minor merging only, wheras non-FGs on 
average form later. The early assembly of FGs leaves sufficient time for galaxies of L ~ to merge into the 
central one by dynamical friction, resulting in the large magnitude gap at z = 0. A fraction of 33±16% of the 
groups simulated are found to be fossil, whereas the observational estimate is ~ 10-20%. The FGs are found 
to be X-ray over-luminous relative to non-FGs of the same optical luminosity, in qualitative agreement with 
observations. Finally, from a dynamical friction analysis is found that only because infall of L ~ L* galaxies 
happens along filaments with small impact parameters do FGs exist at all. 

Subject headings: cosmology: observations - cosmology: - dark matter -galaxies: clusters: general -galaxies: 
formation 



1. INTRODUCTION 

Galaxy groups represent the natural environment where 
galaxy mergers are expected to occur most frequently, be- 
cause of the low velocity dispersion of the group galaxies 
compared to that in galaxy clusters, as well as the large galaxy 
number density. If merging is an efficient process within 
galaxy groups, then there may exist systems in which mem- 
ber galaxies merged to form a large elliptical type galaxy, sur- 
rounded by considerably smaller galaxies, for which the dy- 
namical friction infall time exceeds a Hubble time. In the past 
a lot of work have been done to understand the formation and 
dynamical state and evolution of dense compact groups (Di- 
aferio et al. 1994; 1995; Hernquist et al. 1995; Governato 
et al. 1996). Early numerical simulations of Barnes (1989) 
showed that compact group members can merge to form a 
large elliptical galaxy on timescales considerably less than 
the Hubble time. An observational counterpart to this, called 
'fossil groups' was discovered by the ROSAT X-ray satellite 
(Ponmanetal. 1994). The archetype is RXJ 1340.6+40 18 at 
z=0.171 with a bright isolated elliptical galaxy of magnitude 
Mr ps -23.5 (/i=0.7) surrounded by an extended halo of X- 
ray emitting hot gas, and only much fainter satellite galaxies. 
A fossil group (FG) is defined as a group where a large el- 
liptical galaxy (BG1) embedded in a X-ray halo dominates 
the bright end of the galaxy luminosity function, the second 
brightest group member (BG2) being at least 2 R-band mag- 
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nitudes fainter (Jones et al. 2000, hereafter JPF00). FGs host 
nearly all "field" E galaxies brighter than Mr = -22.5 (/i=0.5) 
(Vikhlinin et al. 1999, hereafter V99) and have mass-to-light 
ratios comparable to those of galaxy clusters, M/L R ~ 300. 
Their X-ray luminosity is for a given X-ray temperature on 
average larger than for normal groups, such that they appear 
to fall on the extention of the cluster Lx - Tx relation to lower 
temperatures (Jones et al. 2003, hereafter J03). With a num- 
ber density of ~ 2xl0~ 7 Mpc" 3 -4xl0~ 6 , they are not rare: they 
constitute ~10-20% of all clusters and groups with an X-ray 
luminosity greater than 2.5 ■ lO 42 /^ ergs" 1 (J03, V99). The 
luminosity function of FGs shows a lack of L* galaxies al- 
though such groups have masses of typically 25% of the mass 
of Virgo cluster. Whereas the Virgo Cluster, e.g., contains 
six L* galaxies, FGs contain typically none, apart from the 
central, bright E galaxy. In terms of the cumulative substruc- 
ture function, which is the number of objects with velocities 
greater than a specific fraction of the parent halo's veloc- 
ity, FGs like RXJ1340.6+4018 show substructure functions 
similar to the Milky Way and M31 rather than, e.g., Virgo 
(D'Onghia & Lake 2004). Hence FGs are excellent probes for 
testing the present understanding of the cosmological struc- 
ture formation, also because X-ray observations strongly in- 
dicate that they are well relaxed, virialized systems. 

In this Letter we present the first attempt, based on self- 
consistent, cosmological simulations, to address the physical 
processes that lead to the formation of FGs. 

2. THE CODE AND NUMERICAL EXPERIMENTS 

We performed simulations of 12 galaxy group-sized dark 
matter haloes in the low-density, flat cold dark matter 
(ACDM scenario) with f2 M =0.3, ft A =0.7, h =H /100 km s" 1 
Mpc"'=0.7 and 08=0.9. The simulations were performed us- 
ing the TreeSPH code briefly described in Sommer-Larsen 
et al. (2005). The code incorporates the "conservative" en- 
tropy equation solving scheme of Springel & Hernquist 2002; 
chemical evolution with non-instantaneous recycling of gas 
and heavy elements tracing 10 elements (H, He, C, N, O, 
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Mg, Si, S, Ca and Fe; Lia, Portinari & Carraro 2002a,b); 
atomic radiative cooling depending on gas metal abundance 
and a redshift dependent, meta-galactic UV field; continuing, 
strong galactic winds driven by star-bursts (SNII), optionally 
enhanced to mimic AGN feedback. A fraction fw of the en- 
ergy released by SNII explosions goes initially into the ISM 
as thermal energy, and gas cooling is locally halted to repro- 
duce the adiabatic super-shell expansion phase; a fraction of 
the supplied energy is subsequently (by the hydro code) con- 
verted into kinetic energy of the resulting expanding super- 
winds and/or shells. 

The original dark matter (DM)-only cosmological simula- 
tion was run with the code FLY (Antonuccio et al. 2003), for 
a cosmological box of 150/r'Mpc box side length. When 
re-simulating with the hydro-code, baryonic particles were 
"added" adopting a global baryon fraction of fb = 0.12. The 
mass resolution was increased by up to 2048 times, and the 
force resolution by up to 13 times. The initial redshift for 
the cosmological run was z,=39. We randomly selected 12 
groups for re-simulation. The only selection criterion was that 
the groups should have virial masses close to (within 10%) 
lxlO 14 M Q , where the virial mass is the mass at z=0 inside the 
virial radius, defined as the region for which the average mass 
density is 337 times the average of the Universe. The corre- 
sponding virial radius is about 1 .2 Mpc, and the virial (X-ray) 
temperature is about 1.5 keV. The purpose of this project was 
to study a cosmologically representative sample of groups, so 
no prior information about merging histories, was used in the 
selection of the 12 groups. Particle numbers were about 2.5- 
3xl0 5 SPH+DM particles at the beginning of the simulations 
increasing to 3-3.5xl0 5 SPH+DM+star particles at the end. A 
novelty was that each star-forming SPH particles of the ini- 
tial mass is gradually turned into a total of 8 star-particles. 
This considerably improves the resolution of the stellar com- 
ponent. SPH particles, which have been formed by recycling 
of star-particles, will have an eigth of the original SPH parti- 
cle mass — if such SPH particles are formed into stars, only 
one star-particle is created. As a result the simulations at the 
end contain star-particles of mass m*=3.1xl0 7 /i _1 M , SPH 
particles of masses m gas =3.1xl0 7 and 2.5x10 s /t'Mq, and 
dark matter particles of ;7?dm=1-8x10 9 /i _i M q . Gravitational 
(spline) softening lengths of 1.2, 1.2, 2.5 and 4.8 /i _1 kpc, re- 
spectively, were adopted. 

To test for numerical resolution effects one of the 12 
groups was in addition simulated at eight times (4 for star- 
particles) higher mass and two times (1.6 for star-particles) 
higher force resolution, than the "standard" simulations, 
yielding star-particle masses m*=7.8xl0 6 h~ l M Q , SPH par- 
ticle masses m gas =7.8xl0 6 and 3.1xl0 7 /i _I M , dark mat- 
ter particle masses todm=2.3x10 8 h~ l M Q , and gravitational 
(spline) softening lengths of 0.76, 0.76, 1.2 and2.4/r'kpc, re- 
spectively. For this simulation particle numbers are ~ 1 .4xl0 6 
SPH+DM particles at the beginning of the simulation increas- 
ing to 1.6xl0 6 SPH+DM+star particles at the end. In pre- 
vious simulations of galaxy clusters (Romeo et al. 2005a,b, 
Sommer-Larsen et al. 2005) we have found that in order to 
get a sufficiently high ICM abundance a combination of a 
large value of fw and a fairly top-heavy initial mass func- 
tion (IMF) has to be employed. We used for the present 
simulations: f w =0.S, an IMF R <^|he r Arimoto-Yoshii type. 

For each of the twelve simulated galaxy groups we com- 
puted the luminosity function of the group member galaxies 
in the B, R, J and K bands (Romeo et al. 2005a). Our group 
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FIG. 1. — The composite group LF for the fossil groups (triangles) and 
non fossil groups (pentagons) of our sample. The observational LF for 
RXJ1340.6+4018 (stars) is shown for comaprison (JPF00). 



sample is divided into two classes according to the magni- 
tude of the R band gap between the brightest galaxy and 
the second brightest galaxy: FGs (Am^s > 2 and non-FGs 
(Am 12 . R < 2). 

A well known problem in the modelling of galaxy groups 
and clusters, is the development of late-time cooling flows 
with bases at the position of BG1 and associated, central star- 
formation rates, which are too large compared to observa- 
tions. In calculating the optical properties of the group galax- 
ies we correct in a crude way for this by removing all stars 
formed at the base of the cooling flow since redshifts z con = 2 
or 1 . Both redshifts correspond to times well after the bulk of 
the group stars have formed. The correction amounts to 20- 
40% in terms of numbers of BG1 stars. Using z colT = 2 or 1 
results in minor differences, so we adopt z C on-=2 in this paper 
— for further discussion of this point see Sommer-Larsen et 
al. (2005). In cases where a similar correction of BG2 is ap- 
propriate (typically for non fossil groups, where BG2 enters 
into the main dark matter halo fairly late, z ^0.5) such a cor- 
rection is applied to BG2. These corrections are quite minor, 
<10% in terms of numbers of stars. 

Figure 1 shows the composite group luminosity function 
(LF) for the 4 FGs (triangles) and 8 non-FGs (pentagons) of 
our sample, adopting a bin size of one magnitude. The large 
magnitude gap Am12.fi between the brightest galaxy and the 
second brightest galaxy is clear from this figure for the FGs 
compared to the absense of a gap for the non-FGs. The ob- 
servational LF (shown with stars) is taken from JPF00. The 
shape of the LF at the bright end, and the typical number of 
member galaxies of our FGs is consistent with the observa- 
tions to our resolution limit of Mr ~ -18.5 (see below). Note 
though, that the LF of RXJ 1340.6+40 18 (JPF00) is only pho- 
tometric, so a LF with members spectroscopically confirmed 
(Mendez de Oliveira et al. 2005) would enable a more fair 
comparison. In order to test how resolution affects our results 
we carried out one comparison simulation at higher resolution 
(Sec. 2). The agreement between the luminosity functions of 
the "standard" and higher resolution run is excellent above the 
resolution limit of the normal resolution run, Mr ~ -18.5, as 
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FIG. 2. — Ami2,R of each simulated group as a function of formation time 
as defined as the epoch in which the group assembled 50% of the system's 
final mass 



will be discussed in D'Onghia et al. (2005, paper II - see also 
Romeo et al. 2005a for details). 

Figure 2 displays an interesting correlation between forma- 
tion time of the group, defined as the epoch at which 50% of 
the system's final virial mass is assembled, and the magnitude 
gap Am\2.R . This correlation indicates that FGs are systems 
formed at z> 1 , and evolved relatively quietly until the present 
epoch, whereas non-fossil groups are generally late forming. 
The earlier a galaxy group is assembled, the larger is the mag- 
nitude gap in the R band, Am 12jR , at z=0. 

A characteristic feature of observed FGs is an excess of 
X-ray luminosity with respect to "normal" galaxy groups 
of the same optical luminosity. This feature is reproduced 
in our simulations: the FGs with mean optical luminosity 
L s = (1.28±O.O6)xlO n L haveL x = (6.3±l.l)xl0 43 erg-.T 1 
whereas the non-FGs with L R = (1.26±0.04)xlO n L Q have 
median X-ray luminosity Lx = (1 .7 ± 0.6)xl0 43 erg-i" 1 show- 
ing a clear offset between the typical X-ray luminosity of FGs 
compared to non-FGs, in qualitative agreement with observa- 
tions. The reason for this is likely that as the simulated FGs 
are systems which form early, gas in the inner regions has 
more time to cool, relative to non-FGs, resulting in stronger 
central cooling flows, and hence higher X-ray luminosities. 
Note, however, that stronger than observed cooling flows at 
late times occur in general at the centers of the simulated 
groups. This results in low central gas entropies and X-ray 
luminosities about a factor of 3 larger than observed. Similar 
offsets are found for the L^-Tx relation of simulated galaxy 
clusters, when an IMF of the Arimoto-Yoshii type is adopted 
(Romeo et al. 2005b). 

4. DISCUSSION 

We have found a correlation between the magnitude gap be- 
tween the BG1 and BG2 and the formation time of the dark 
halo of the groups considered. What is the origin of the large 
magnitude gap at the bright end of the LF of FGs? In Figure 
3 we show the evolution with redshift of the stellar mass of 
BG1 and the second brightest galaxy at any given time, BG2, 
for a FG and a non-FG from our sample. At late times (z <1) 
the "magnitude" gap in the FG increases due to minor merg- 
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FIG. 3. — The evolution of the stellar mass of BG1 (filled circles) and BG2 
at any given time (filled triangles) for a FG. Note that for the FG the BG2 
is the second brightest object at each output time. The BG2 merges into the 
BG1 at z=2.2 and the last major merger of a BG2 into BG1 occurs at z=1.4 
(indicated by arrows). For a non-FG is shown M*(z) of the galaxy which 
becomes the final BG1 (BG2) by open circles (triangles). 



ing/stripping. We interpret this result as follows: since the 
dark halo is assembled at z > 1, there is plenty of time for the 
BG2s to first be tidally stripped during close passages to the 
BG1 and finally merge into the BG1, by dynamical friction. 
Note that in Figure 3 for the FG the BG2 is the second bright- 
est object at each output time. At z=2.2 the BG2 at that time 
merges into the BG1. The last major merger of the BG2 into 
BG1 occurs at z= 1.4 (see arrows in Figure 3). For the non-FG, 
the stellar mass of both BG1 (empty circles in Figure 3) and 
BG2 (empty triangles) mostly increases with time till z=0. In 
this case the final halo assembles its mass fairly recently, and 
the two largest, more similar sized galaxies have not had time 
to merge. 

We now analyze whether the fossil group phenomenon can 
be explained by dynamical friction effects. First, we compute 
the time t- m f it takes for a satellite galaxy moving in a main 
dark matter halo to be accreted into the central region of the 
halo: we assume singular isothermal sphere models for both 
the main halo and the satellite, with characteristic circular ve- 
locities of V# and V$. If we assume that the satellite initially 
moves in a circular orbit (one can show that this gives the min- 
imum infall timescale for a given initial impact parameter) at 
a speed V = Vh and initial radius ro, and that f,-„y is consid- 
erably larger than the dynamical timescale, the satellite will 
feel a small torque which slowly reduces its orbital angular 
momentum, keeping the orbit almost circular while it spirals 
towards the centre of the halo: 



dl 



dr 



dt dt 



(1) 



denoting / the specific angular momentum of the satellite, r 
its orbital radius and aj„, the drag acceleration opposed to 
the direction of motion. The latter can be expressed by (e.g., 
Binney & Tremaine 1987): 

47rln AG 2 mp 
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(2) 



where X = V/\/2, <th = V /Vh, m is the mass of the satellite, p 
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the (mainly dark) matter density, and In A the Coulomb loga- 
rithm, which from numerical simulations is found to be <~ 1 .5. 

The radius of the satellite halo is set by tidal truncation r s = 
{mj'hM) x ^r, where M is the mass of the main halo inside of 
r. For a isothermal sphere model we have r s = -^-r and 

m = ^ - r, inserting these equations and eq.(2) in eq.(l) we 
obtain: 

^-|- i2 - 4 100kpcV700km/sJ V250km/sJ ^ ' 

(3) 

for In A=1.5. Hence, an L» satellite (Vs ~ 250km/s) being de- 
posited as close as 100 kpc from the BG1 (given that the group 
virial radius at z=0 is 1.2 Mpc) within a group of Vr ~ 700 
km/s (corresponding to a virial mass of 10 14 M Q ) on a low- 
eccentricity orbit has an infall time comparable to the Hubble 
time. It then at first seems puzzling that fossil groups ex- 
ist at all. However, the infall time is proportional to ro, and 
as will be discussed in paper II in the simulations infall hap- 
pens along filaments with initial impact parameters as small 
as 5-10 kpc — this makes infall times of L* galaxies less than 
the Hubble time, and is the main reason why fossil groups 
exist. If the early merging history is responsable of the for- 
mation of fossil groups, why are there no clusters with a sim- 
ilar lack of L ~ L* galaxies as a result of cosmic variance? 
As characteristic initial impact parameters scale linearly with 
Vh, fixing Vg we obtain oc V^. For a typical rich clus- 
ter V H ~ 1400-2100 km/s, so it follows that the typical in- 
fall time for a L* galaxy will be 10-30 larger than in groups. 
Hence massive galaxies have had no time to merge with the 
central object yet. 

Given that merging is such an important ingredient in FG 
formation, we check that merger rates are correctly estimated. 
When calculating the optical luminosities of the BGls (and in 
some cases also BG2s) we correct for the effect of late time 
cooling flows by not including stars formed later than z C orr=2 
at the base of the cooling flow. These stars will, however, still 
contribute to the gravitational field of BG1, which might lead 
to an increase of the merger rate. To test for this we ran one 
simulation, where star particles formed at the base of the cool- 
ing flow (at BG1) were removed entirely from the simulation. 
This did not lead to any significant change of the merger rates. 
Moreover, as discussed above, analysis of our high-resolution 



run indicates that our simulations are not affected by "over- 
merging". 

Our work shows that in the hierarchical structure formation 
scenario FG halos assemble earlier then non-FG halos, and 
they preferentially reside fairly isolated in low-density envi- 
ronments. Once the halos assembled 50% of their final mass, 
they evolve quiesciently. Consequently, they experience only 
little late infall of galaxies that could stop merging and fill the 
magnitude gap of the LF again. 

We now consider the regions around the two FGs with the 
largest magnitude gaps, Am12.fi > 3, compared to the regions 
around the non-FGs: between 2 and 5 Mpc, the region around 
FGs shows a total matter density in units of the mean of 
the Universe at z=0, p tot /p = 2.06 ±0.23, whereas the region 
around non-FG halos has p t0 , jp = 4.86 ± 0.99. So the region 
within 5 Mpc around FGs is found to be of lower density (at 
about the 3-a level), compared to non-FGs. For the region 
2<R<7.5 Mpc, the corresponding numbers are 1 .88 ± 0.39 
and 3.11 ±0.45, so the region around the FGs is again under- 
dense with respect to that of the non-FGs (at about the 2-a 
level). From R^IO Mpc and outwards there is no statistically 
significant difference between the environments and the two 
types of groups. This prediction can be observationally tested. 
If FGs are X-ray emitting systems in place already at z=l, 
they might be observed with similar techniques used to dis- 
cover high-redshift galaxy-clusters (Mullis et al. 2005). Once 
the extended X-ray source is identified in the XMM-Newton 
archive, subequent VLT multi-object spectroscopy combined 
with optical/NIR imaging could reveal the possible presence 
of FGs at high redshift. Furthermore, the correlation function 
of X-ray emitting galaxy groups cross-correlated with galax- 
ies in the nearby universe, e.g. from SDSS data, could be 
used to test whether FGs preferentially populate low-density 
regions of the Universe. 

Of our 12 simulated groups, 4 become FGs corresponding 
to a fraction of 33±16%, larger but still comparable to the ob- 
servational estimate of 10-20%, although the sample of simu- 
lations and observational data is statistically limited. 

We are grateful to the referee for valuable suggestions and 
to F. van den Bosch, G. Lake, D. Pierini and TPonman for 
discussions. 
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